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During spontaneous imbibition a wetting liquid Is drawn Into a porous 
medium by capillary forces. In systems with comparable pore length 
and diameter, such as paper and sand, the front of the propagat- 
ing liquid forms a continuous Interface. Sections of this Interface 
advance in a highly correlated manner due to an effective surface 
tension, which restricts front broadening. Here we investigate water 
imbibition in a nanoporous glass (Vycor) In which the pores are much 
longer than they are wide. In this case, no continuous liquid-vapor 
interface with coalesced menisci can form. Anomalously fast imbibi- 
tion front roughening is experimentally observed by neutron imaging. 
We propose a theoretical pore network model, whose structural de- 
tails are adapted to the microscopic pore structure of Vycor glass, 
and show that it displays the same large scale roughening charac- 
teristics as observed in the experiment. The model predicts that 
menisci movements are uncorrelated. This indicates that despite 
the connectivity of the network the smoothening effect of surface 
tension on the imbibition front roughening is negligible. These re- 
sults suggest a new universality class of imbibition behavior which 
Is expected to occur in any matrix with elongated. Interconnected 
pores of random radii. 

liquid imbibition | porous media | interface roughening | neutron imaging | 
pore network model | computer simulations 

Abbreviations: NVG , nanoporous Vycor glass 

Many everyday processes involve the flow of a liquid into 
a porous matrix, for instance when we dunk a biscuit 
into coffee, clean the floor with a cloth, or get drenched with 
rain. The same process is also important in nature, e.g. for 
water to reach the tips of the tallest trees or to flow through 
soil, and crucial for different industrial processes, ranging from 
oil recovery and chromatography to food processing, agricul- 
ture, heterogeneous catalysis, and impregnation (for reviews 
see PEmil). 

The above processes are examples of imbibition (Fig. [T]). 
Imbibition of a liquid into a porous matrix is governed by 
the interplay of capillary pressure, viscous drag, volume con- 
servation and gravity. The porous matrix often has a com- 
plex topology. The inhomogeneities result in variations in the 
local bulk hydraulic permeability and in the capillary pres- 
sure at the moving interface. Nevertheless, the invasion front 
during solely capillarity-driven (i.e. spontaneous) imbibition 
advances in a simple square-root-of-time manner, according 
to the Lucas- Washburn law [5] [6]. This is a result of the 
time-independent mean capillary pressure and the increasing 
viscous drag in the liquid column behind the advancing front. 
It is valid down to nanoscopic pore sizes [T] [51 [9] and par- 
ticularly robust with regard to the geometrical complexity of 
the porous matrix [I] O IIOI ill] . The evolution of the inva- 
sion front displays universal scaling features on large length 
and time scales, which are independent of the microscopic de- 
tails of the fluid and matrix [El [111 [H [I5l [H [iZl [18] . This 
parallels the elegance of critical phenomena. 



Typically imbibition is studied using paper [141 1151 116] 
or Hele-Shaw cells [171 118| . In these systems, pore space is 
laterally highly interconnected. This results in a continuous 
liquid-gas interface, whose advancement is spatially correlated 
due to an effective surface tension [19| . Consequently menisci 
advancement beyond the average front position is slowed down 
whilst menisci lagging behind are drawn forward. Hence, the 
roughening of the interface is slowed down. By contrast, in 
many real porous systems, e.g. rock and soil |20) . the pore 
network consists of elongated pores with reduced connectivity 

Here we investigate the spontaneous imbibition of water 
into nanoporous Vycor glass (NVG), which is a silica substrate 
with an interconnected network of nanometer-sized, elongated 
pores |21) . The narrow pores lead to capillary pressures of 
several hundred times atmospheric pressure. This means that 
gravity would only halt capillary rise after several kilometers 
and several billion years |22| . Hence, with this system, we are 
able to observe pure spontaneous imbibition over large length 
(centimeter) and long time (hours) scales. The observation 
of the advancing front is difficult since it is deeply buried in- 
side the matrix |23l 124] . Nevertheless, neutron radiography 
|25U26lf27] allows us to image the liquid inside porous materi- 
als |28l 129] . Recent technical improvements in neutron imag- 
ing provide the spatial and temporal resolution (tens of fim 
and tens of s) to follow imbibition and to obtain quantitative 
information on the morphological evolution of the progressing 
interface. 

Our experiments and simulations show that the interface 
roughness as measured by the interface width w{t) increases 
much faster than observed previously, namely w{t) oc with 
P ~ 0.45. This is close to the square-root-of-time progres- 
sion of the invasion front H(t) ~ t^'^ . The propagation front 
hence comprises an almost constant fraction, w{t)/H{t), of 
the occupied part of the matrix, including voids and over- 
hangs. We find that lateral correlations of the invasion front 
are short-ranged and independent of time. This means that, 
in the present case, surface tension is irrelevant in a coarse 
grained description of interface broadening on macroscopic 
length scales. Similar broadening has been observed, but for 
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difFerent experimental conditions, namely drainage and forced 
imbibition of less- wetting fluids |30l 131] . 

Results and Discussion 

Quantitative Experimental Characterization of the 
Imbibition Front. We investigate imbibition in NVG which 
contains pores with a mean radius rav ~ 4 nm, a radius poly- 
dispersity of 20%, a pore aspect ratio a = L/2rav between 5 
and 7 (where L is the pore length), and a porosity of about 
30% 13 [211132] [33]. (For details see Methods.) The bottom 
face of an empty NVG is brought into contact with the surface 
of a water reservoir. Capillary forces draw the liquid into the 
porous matrix. 

The invasion front appears as a bright region when imaged 
using reflected light (Fig. [2] top row). The pronounced (mul- 
tiple) light scattering is caused by filled and empty regions 
randomly alternating on the length scale of the wavelength 
of visible light; several hundred nanometers |34) . These spa- 
tial scales are much larger than any intrinsic length of the 
porous matrix, such as the pore diameter or the pore-pore 
distance. Thus strong scattering of light at the advancing in- 
terface indicates structures of filled and empty parts (voids) 
of surprisingly large size. It prevents also a quantitative de- 
termination of the liquid profile in the propagating front and 
hence of the width of the interface. By contrast the spatial 
and temporal resolutions of neutron imaging [35| allow quanti- 
tative measurements in systems such as NVG (Fig. [2] bottom 
row). 

From the neutron images we determine the spatial and 
temporal evolution of the local filling degree < f{x, z, t) < 1. 
Due to the projection in the y-direction, this is the average 
amount of filled pore space at lateral position x, height z 
and time t. Its lateral average, that is the vertical concentra- 
tion profile f{z,t) = {f{x,z,t))x, is shown in Fig. [3^. The 
time-dependence of the front height, quantified by the mean 
median rise level H{t) = (5;(/=0.5, x, t))x, follows the Lucas- 
Washburn \/t-law (Figs. [3K,b, solid lines), consistent with 
previous studies Fits of Gauss error functions to the pro- 
files yield the time-dependence of the width w{t) (Fig. 
The fit of w{t) oc results in a growth exponent of the width 
or roughness, /3 = 0.46 ± 0.01 (Fig. [3J:, solid line). The value 
/3 = 0.46 significantly exceeds previous theoretical predictions, 
in particular those from phase-field models which are based on 
quenched, random fields. Such models predict slower rough- 
ening dynamics with /3 ~ 0.19 and a strong spatial correlation 
of the height fiuctuations within the moving interface [13| . 

Instead of the median rise level averaged in x direc- 
tion, H{t) — {z{f —0.5, x,t))a:, we now consider the local x- 
dependent median rise level h{x,t) = z{f =0.5, x,t) (Fig. [SJd) 
to investigate fluctuations in x-direction, i.e. within the front. 
We calculate the height-height correlation function: 

C{£,t) = ^/{{h(x,t) - h{x+£,t))^)^ . [1] 

The observed fluctuations in C{£, t) (Fig. [3jl) are due mainly 
to the limited data density and stray gamma radiation from 
the reactor and instrument hitting the camera. The data 
exhibit neither scaling of C{l,t) with I nor any indication 
of spatial correlations in the experimentally accessible range 
75 Jim < £ < 4000 pm. Although the correlations are reduced 
due to the projection in y-direction, the absence of any de- 
tectable correlation is in contrast to all previously reported 
experiments and theories on imbibition front roughening. 

Pore-network Model. 

No theoretical model is available that is consistent with 
our system and which predicts the spontaneous imbibition be- 



havior observed. An ensemble of independent pores of random 
but constant radius exhibits a roughening exponent 1/2 since 
the meniscus heights evolve independently from one another 
as hi{t) = ai\/t with random pre-factors ai. However, this in- 
dependent pore model is inappropriate for Vycor glass, since 
the pore radii vary strongly along individual pores (see Fig. 
1). An ensemble of independent pores with radii which vary 
randomly along their length has a roughening exponent of 1/4 
(see Appendix) , which does not agree with the experimentally 
observed value. Thus independent pore models do not explain 
the observed exponent. A roughening exponent /3 « 1/2 has 
recently been reported within the framework of a lattice gas 
model for spontaneous imbibition |41| . This model is appro- 
priate for silica aerogels with an extremely large porosity of 
87-95 % and gives rise to a continuous liquid-gas interface. 
Consequently one expects here an effective surface tension to 
be present, inducing height-height correlations in the advanc- 
ing imbibition front. The model details are thus not appropri- 
ate for NVB. To our knowledge, all other existing theoretical 
models (for an overview see U) are also incompatible with 
our experimental observations which are i) fast broadening 
dynamics with a growth exponent close to 1/2, ii) absence of 
height-height correlations in the advancing imbibition front. 

Hence we propose a pore-network model |36l I37| adapted 
to our experimental situation. This consists of individual, 
elongated capillaries arranged in a two-dimensional square 
lattice with laterally periodic boundary conditions. Capil- 
laries are connected at nodes and inclined at 45°. All cap- 
illaries have the same length L, while the radius r of each 
capillary is randomly chosen from a uniform distribution with 
mean radius rav and width 25^, i.e. disorder strength Si/ra_v 
|21l 133). (For details see Methods.) We investigated aspect 
ratios 2.5 < a < 10 and polydispersities 0.1 < (5r/fav < 0.4. 

The pressure at all nodes at the bottom of the lattice are 
set to zero while at the menisci the Laplace pressure prevails. 
This pressure difference drives the flow through the capillar- 
ies. This flow is opposed by viscous drag according to Hagen- 
Poiseuille's law. During the whole process, volume conserva- 
tion must be maintained. When a meniscus reaches an empty 
node it 'jumps' over the node, generating new menisci in a 
distance S = L/lOO from the node (Fig. [4^). This implemen- 
tation of node crossing avoids a microscopic treatment of the 
flUing process of the nodes and is valid as long as this is not the 
rate limiting step, which we will discuss below. If the Laplace 
pressure of a meniscus exceeds the node pressure the meniscus 
is arrested at the distance S until the node pressure increases 
beyond the Laplace pressure (Fig. \^). Then propagation of 
the meniscus resumes. 



Computer Simulation Results of Imbibition. With 
our model, we observe a strong roughening of the imbibi- 
tion front (Fig. [4]1) with fast moving menisci advancing 
through sequences of thin capillaries and arrested menisci 
lagging behind. Quantitatively, the computer simulations 
yield a mean rise level H{t) = {hi{t)} (where {. . .) denotes 
an average over all menisci labelled by the index i), which 
obeys the Lucas- Washburn ^/t-behavior (Fig. \5^). The width 

^(^) ~ ii^^i't)) ~ i^ii^))^)^^^ increases rapidly as w{t) oc 
with /? = 0.42 ±0.01 or 0.45 ±0.01 for the smallest and largest 
polydispersities, respectively (Fig. [SJj). We find a slight up- 
ward trend of 'w{t) at large times (Fig. [5J3, inset), which is 
also suggested in the experimental data (Fig. [3]:) and indi- 
cates that the asymptotic value of the growth exponent might 
be larger. This is consistent with the fact that a smaller /? 
is observed when the asymptotic behavior is approached later 
(as in the case of the smaller polydispersity). This implies that 
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the asymptotic value is closer to that found for the larger poly- 
dispersity with an increased uncertainty, i.e. /3 = 0.45 ± 0.02. 
Variation of the aspect ratio in the range 2.5 < a < 10 gave 
identical results for the exponent /3. 

We systematically studied finite size effects, especially on 
Remarkably, we only find a dependence on the lateral 
system size A'x for the smallest system size A'x = 4 (Fig. [5j;). 
This provides an upper bound for the characteristic length 
scale ^{t) of height-height correlations. Within the framework 
of the scaling theory of roughening [381139) . the interface width 
in a finite system of lateral size Nx is expected to behave as 



WN^ (t) 



const. 



for 
for 



m/L « Nx 
m/L » Nx 



|2] 



The data (Fig. [5t) suggest ^(i) < 4L, implying that the 
roughening dynamics are not or are only weakly spatially cor- 
related. This is confirmed by the height-height correlation 
function C{£,t) (eq. [T] Fig. HJl), which saturates quickly 
(around £ m L). Scaling theory [381 139] predicts saturation 
of C{£,t) for £ > ^{t), which implies ^{t)/L = 0(1) indepen- 
dent of time t. This finding supports our experimental result 
(Fig. [3ji) that any spatial roughness correlations are absent 
and extends its validity down to pore-pore distances and thus 
towards the nanometer-scale, i.e. far beyond our experimental 
resolution. 

Experiments and Simulations Both Yield an 
Anomalous Roughness Growth Exponent. Experiments 
and simulations exhibit corresponding behavior, even on a 
quantitative level, i.e. progression of the imbibition front ac- 
cording to the Lucas- Washburn y^-law, fast broadening of the 
front with a large growth exponent j3 « 0.45 and short range 
height-height correlations over a maximum of 1-2 pore lengths 
only. The pore-network model with its elongated pores hence 
successfully mimics the characteristics of NVG. In such mor- 
phologies, all menisci are restricted to individual pores and 
thus cannot interact via an effective surface tension. In the 
absence of interactions between individual menisci local pro- 
cesses at the junction become important for the front rough- 
ening. The dynamics of junction filling, as analyzed in [40] . 
describes a threshold mode in which meniscus propagation is 
halted whilst the junction is filled and the new menisci in the 
adjacent pores form. This takes only a few milliseconds in 
nanometer sized pores with filling heights up to 2 cm (as in 
the present case), for which gravity is negligible. Once the new 
menisci have formed, those in the thicker pores are arrested 
as long as their Laplace pressure is larger than the pressure 
within the junction. These arrests can last much longer than 
the filling process, up to times of the order of the age of the 
propagation front, which can be several hours in our experi- 
ments. Thus for the asymptotic (long time) behavior of the 
broadening dynamics, the filling process of individual junc- 
tions is negligible and front roughening is mainly influenced 
by the arrests after the filling of the junctions. As a conse- 
quence, the distribution of pore diameters and the frequen- 
cies of junctions is expected to be much more important than 
the topology, in particular the dimensionality, of the network. 
This contrasts with the role the latter usually plays in surface 
roughening and critical phenomena [381 139[ . 

It should be noted that it is crucial that the pores in NVG 
are connected [21U33| . since an ensemble of independent pores 
with randomly varying radii over their lengths would give a 
roughness exponent 1/4 (see Appendix). Perhaps counterin- 
tuitively, the introduction of junctions, i.e. branch points or 
crossings, enhances height fluctuations and hence the front 
width. This is because at each branch point one meniscus can 



split into two (or more) upwards moving menisci with one typ- 
ically moving faster than the other, or even one moving and 
the other stopping until the node pressure exceeds its Laplace 
pressure. 



Conclusions 

Pore Aspect Ratio Determines the Universality Class. 

Our results show that spontaneous imbibition crucially de- 
pends on the pore aspect ratio a. For short pores (small a), 
neighboring menisci coalesce and form a continuous imbibi- 
tion front. Thus the smoothening effect of an effective surface 
tension within the interface leads to a slow broadening of the 
front. Various theoretical models describe the roughening of 
a vapor-liquid interface during spontaneous imbibition in the 
presence of an effective surface tension, e.g. phase fleld mod- 
els [19| . These predict a roughening exponent /3 ~ 0.19. The 
elongated pores in nanoporous Vycor glass (large aspect ratio 
a) inhibit the formation of a connected vapor-liquid interface. 
In this case the individual menisci cannot interact via an effec- 
tive surface tension and the broadening of the imbibition front 
is anomalously fast with /3 « 1/2 establishing another univer- 
sality class. The regime of weak roughening (small a) must be 
separated from the regime of strong roughening (large a) by a 
critical value Uc of the aspect ratio, but its precise value will 
depend on structural details of the pore network, in particular 
the pore junction geometry. 

We want to stress that strong imbibition front broadening 
is not linked to the nanometer size of the pores. However, its 
experimental observation over large length and time scales sig- 
nificantly benefits from the dominance of capillary forces over 
gravitational forces, which results from the nanometer-sized 
pores. The theoretical model employs macroscopic hydrody- 
namic concepts only. Therefore, strong interfacial broaden- 
ing is a consequence of any spontaneous imbibition process 
in porous structures with interconnected elongated capillaries 
independent of their macroscopic extension and mean pore di- 
ameter. It is not only important to nanofiuidics, but for liquid 
transport in porous media in general. 

Our observation of a new universality class of strong in- 
terfacial broadening is thus a very general finding, which has 
been made possible due to recent improvements in the reso- 
lution of neutron imaging [35| . The front roughness is crucial 
for many processes, such as water transport in geology, flux 
in oil recovery, glueing, dying and impregnation. Our results 
enable us to link the broadening dynamics during these pro- 
cesses to the properties of the porous materials. To what 
extent this behavior can be described with alternative models 
for transport in porous media, e.g. models which consider a 
saturation-dependent hydraulic permeability of the pores [4], 
warrants further investigation. 

Materials and Methods 



Neutron Imaging. 

The nanoporous Vycor glass (NVG) consists of an interconnected networl< 
of elongated pores witii a mean radius Tav = 4 nm, a radius polydispersity 
Sr/fav = 0.2, a pore aspect ratio 5 ^ a ^ 7, and a porosity of about 30% 
[l||2l][32]|33]. Tiie macroscopic dimensions of the sample are 4.6 X 4.6 X 20 mm^ . 
Its faces, except the bottom face, are sealed to preclude liquid evaporation. To initiate 
imbibition, the bottom face of the sample is brought into contact with the surface 
of a water reservoir. During imbibition, the huge capillary pressure highly compresses 
entrapped air v*/hich is subsequently dissolved in w/ater and hence does not affect our 
experiments. All experiments are performed at room temperature. 
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The neutron imaging experiments are performed at the ANTARES beamline of 
the research reactor FRM II of the Technical University Muenchen (Garching, Ger- 
many) |42| . A beam of cold neutrons passes through an aperture with size D and, 
after a distance L, 'illuminates' the sample, which is situated d — 30 mm in front 
of the scintillator. The geometrical resolution is dli^L/D^ = 75 /im. The trans- 
mitted neutrons are detected using a very thin 'Gadox' scintillator, which does not 
limit the geometrical resolution, and a CCD camera with pixel size 15.97pm. Series 
of images are recorded for total measurement times up to several hours. Individual 
measurement times are 30 s and data transfer times 10 s. In the first few kinetic 
images smearing occurs due to the front moving a significant distance during the 
individual measurement times. However, after about 1000 s, the smearing due to the 
limited time resolution is negligible compared to the spatial resolution (and only these 
data are used for fitting w{t)). Raw images were corrected for detection efficiency, 
background and noise, while corrections for scattered neutrons are not necessary |43| . 

The experimentally determined neutron transmission T{x, Z, t), that is the ra- 
tio of transmitted intensity /(x, 2, and incident intensity Iq{x^ Z^t^ , is related 
to the absorption coefficient S(^X^ 2, t) by 



I{x,z,t) 



-S{x,z,t) d 



[3] 



h){x,z,t) 

where d = 4.6 mm is the sample thickness. The absorption coefficient 

S{x,Z,t) = Sn-iix,Z,t) -\~ f{x,Z,t) Sv, [4] 

depends on the absorption coefficient of the porous matrix Sn-i{x, Z^t^, experi- 
mentally determined from the dry matrix, and on that of the liquid i^wp determined 
from the completely filled matrix providing Sm(^X, Z^t^-\-Sw The filling factor 
f{x, Zf t) can then be determined from the experimentally determined transmission 
T\X^ Z^ t). While silica, and thus NVG, is almost transparent to neutrons, the neu- 
tron beam is strongly attenuated by hydrogen in the water. The contrast is further 
enhanced by the characteristic wavelength distribution of the ANTARES beamline, 
which contains a large fraction of cold neutrons. 

Computer Simulations. The pore-network model consists of capillaries ar- 
ranged on a two-dimensional square lattice inclined at 45 . The system consists of 
A^x and TVz nodes in the horizontal and vertical directions, respectively, with peri- 
odic boundary conditions in horizontal direction. At the nodes, four capillaries are 
connected to each other (Fig. [4}. All capillaries have the same length L, while the 
radius of each capillary is chosen randomly from a uniform distribution with mean 
radius Tav and distribution width (5r, i.e. disorder strength Sr/v^v- We performed 
computer simulations for different lateral system sizes 4 ^ A^x ^ 32 and a vertical 
size up to A^z = 1000, which implies a maximum height H = -\/2 L Nz which 
was not reached by the invasion front within the simulation time. 

The water rises spontaneously from the bottom to the top of the lattice. The 
dynamics are controlled by capillary pressure, viscous drag and volume conservation. 
At each meniscus, i.e. for each capillary J connected to node Z, we calculate the 
capillary pressure given by the Laplace pressure 



c.i 7 

ri 



[5] 



where is the radius of the capillary and cr the surface tension {(J = 72 mN/n- 
for water). Flow through the capillary is driven by the pressure difference A_P^^ = 
Pi — j, where Pi is the pressure at node i. 

According to Hagen-Poiseuille's law, the volume flux Q'^ from node i into cap- 
illary j is 



877 



[6] 



where /l;^ is the length of the liquid column in capillary j of node i and Tj the vis- 
cosity of the liquid [T] = 0.88 mPa s for water). The volume flux determines 
the change of the liquid volume 'V^ and thus of the length of the liquid column 
hi according to = dV^ /dt = n (r^)^ dh^ / dt . Hence, once the node 

pressures Pi are known, the time dependencies of the heights h"^^ (i) are given by 
ordinary differential equations. 

The node pressures Pi are determined by the boundary conditions and volume 
conservation. The boundary conditions are the Laplace pressure at the menisci, P^ ^ , 
and zero pressure at all nodes at the bottom of the lattice which are connected to 



the water reservoir. The volume conservation at each node is given by 



which corresponds to Kirchhoff's law. The sum runs over all capillaries J attached 
to node Z. The resulting set of sparse linear equations is numerically solved to obtain 
the node pressures Pi for a given meniscus height configuration h"^- . The differential 
equations for /l^ are then numerically integrated using an implicit Euler scheme for 
time-stepping. Note that, due to the nanometer-sized capillaries, capillary pressure 
dominates gravity, which can thus be neglected. 

The time step At in the numerical integration of the equations of motion of the 
menisci heights is chosen such that each meniscus moves at most a distance L/IO 
and no meniscus crosses a node. If this would occur for one meniscus. At is reduced 
such that this meniscus reaches the next node and then 'jumps' over the node, gener- 
ating new menisci in a distance S = _L/100 from the node (Fig. [4j), and all other 
menisci are also processed with the reduced At. Similarly, if the meniscus retracts 
due to a negative pressure difference, ^Pl <^ 0, the meniscus is arrested when it 
has approached the node up to a distance 5. Thus a liquid column with a length of 
at least 5 is kept in the capillary, i.e. /i^ ^ (5 always holds. The meniscus is released 
when A_Pj^ > (Fig. [4})). When two menisci meet, they merge and the capillary 
thus is com pletely filled (Fig. which mimics the absence of entrapped air in our 
experimental system. 

During a computer simulation of the time evolution of the model, the aver- 
age rise level H {t^ = {JliitYji of the invading front and its width 'w{t') = 

{{h^{t))i — i^iii))!)^^^ calculated at different times t. Since the invasion 
front contains overhangs and voids, the average (. . .) { is taken over all menisci in- 
dexed by 1. The presented data are averaged over 100 simulation runs using different 
disorder realizations. The statistical error of this average is represented by the error 
bars of the simulation results. 



Appendix: Imbibition front broadening in an inhomoge- 
neous porous medium of independent pores 

For comparison with the proposed model, we here con- 
sider spontaneous imbibition in an ensemble of independent, 
i.e. non-connected or isolated, pores. The radius of a sin- 
gle pore varies randomly with height h such that an ap- 
propriate model for the meniscus motion in such a pore is 
dh/dt = K{h)/h, where K{t) is uncorrelated white noise with 
mean c and variance a, i.e. {ii{h)) = c, — 
a 5{h — h'). For the time T to reach some height H one 
thus gets T{H) = dhh^{h) where ^(h) is white noise with 
mean 1/c and variance a'. Averaging the stochastic variable 
T yields (T) oc (which is the Lucas- Washburn law) and 
for the variance AT^ = ((T - {T)f) oc , which means 
AT oc T^^* . The time to reach height H therefore varies typi- 
cally between + AT and H'^ - AT, vice versa at time T one 
then expects the height H{T) to vary between [T - ATf^ 
and (T -f- AT f^'^ which means AH ^ AT/T^^^ = T^/^. 
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Fig. 1. Schematic representation of spontaneous imbibition of a fluid into a porous matrix. 
The arrows indicate the mean median rise level H{t) and the invasion front width wit). 



Fig. 2. Direct observation of spontaneous imbibition of water into nanoporous Vycor glass 
using visible light (top row) and neutrons (lower row). The reflected light intensity and the local 
liquid concentration, i.e. the filling degree /(x, 2, t), are shown in grayscale and pseudocolors, 
respectively. Snapshots are recorded for 0.1 s (light) and 30 s (neutrons) about every 15 min. 
The lateral direction X and the height z, i.e. the direction of capillary rise, are indicated. The 
width and height of the sample are 4.6 mm and 20 mm, respectively. 
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Fig. 3. Progression and broadening of the imbibition front of water in nanoporous Vy- 
cor glass as quantified by neutron imaging, (a) Laterally-averaged filling degree f[z^t) = 
(/(x, z, £))x as a function of height z and time t. (b) Local median rise level h = 
z{f=0.b,x,t) as a function of lateral position X and time t. (c) Evolution of the front 
width w(t) along with a fit of w oc (solid line). The first few data points show an ap- 
parently increased width due to smearing effects and are thus disregarded in the fit. (See text 
for details.) The inset shows the same data in a log-log representation. Axes units of the inset 
agree with the ones of the main plot, (d) Height-height correlation function C(-^, t) of the 
invasion front at three different times (as indicated). The data are shifted for clarity. In (a) 
and (b) the Lucas-Washburn law ZOCs/t is shown as solid lines. Parts (a) and (b) show only 
about the first half of the collected data, during which most of the front movement occurs. 
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Fig. 4. Computer simulation of spontaneous imbibition in a network of elongated capillaries. 
Rules for the propagation of menisci are designed to mimic the experimental situation, (a) 
After reaching an empty node, the liquid immediately fills the connected capillaries for a distance 
b. (b) After retracting up to (5 toward a filled node, the meniscus is arrested until the 
pressure difference driving the liquid is again positive, (c) When two menisci meet, they merge, 
(d) Snapshots of configurations in a system with aspect ratio a = 5, radius polydispersity 
f^r/Tav = 0.3 and lattice size 16 X 64 at different times (as indicated in units of ns). 
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Fig. 5. Progression and broadening of the imbibition front as observed in computer simu- 
lations based on our pore-network model. Pore aspect ratio a = 5 and standard lattice size 
A^x X A^z = 16 X 1000. Lengths are measured in units of L, times in units of ns. All data 
are averaged over 100 simulations with different disorder realizations, the error bars reflect the 
standard deviation, (a) Mean median rise level = for different polydispersities 

f^r /'"av (as indicated). The lines represent fits of oc \fi.. The inset shows the same data 
in a log-log representation and M oc sft as the dash-dotted line, (b) Evolution of the front 
width for different polydispersities (as indicated). The lines represent fits of oc 

with growth exponents /3 = 0.42 di 0.01, 0.42 di 0.01 and 0.45 di 0.01 for polydispersities 
of 0.1, 0.2 and 0.4, respectively. The inset shows the same data in a log-log representation 
and y/t as the dash-dotted line, (c) Evolution of the front width for a polydispersity 

5rlTa,v = 0.1 and different lateral system sizes A/x (as indicated). The line represents a 
fit of w oc with /3 = 0.42 di 0.01 for A^x > 4. The inset shows the same data in a 
log-log representation, the dashed line represents an upper bound of the width, which is given 
by the difference between the front heights H(t) cc y/t in two homogeneous systems with 
constant minimal and maximal capillary radius r = T^m—^t (slowest front propagation) and 

— '"av+'^r (fastest front propagation), respectively, (d) Height-height correlation C(-^, 
as a function of the distance i at different times t (as indicated) for a sample with polydispersity 
^r/rav = 0.1 and A^x X A^z = 128 X 32. 
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